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Abstract 

We discuss the small- £ behaviour of the next-to-leading BFKL equation, de- 
pending on various smoothing out procedures of the running coupling constant at 
low momenta. While scaling violations (with resummed and calculable anomalous 
dimensions) turn out to be always consistent with the renormalization group, we 
argue that the nature and the location of the so-called hard Pomeron are dependent 
on the smoothing out procedure, and thus really on soft hadronic interactions. 
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Interest in the high energy behaviour of Quantum Chromodinamics [Q has recently 
revived because of the experimental finding |J of rising structure functions at HERA and 
has triggered a number of papers about the next-to-leading (NL) BFKL equation || and 



One of the interesting features of the NL BFKL equation is supposed to be the descrip- 
tion of running coupling effects, which on one hand raises the question of its consistency 
with the renormalization group (R. G.) and on the other hand emphasizes the problem 
of the singular transverse momentum integration around the Landau pole, sometimes 
referred to as the IR renormalon problem |J. 

The purpose of this note, which is based on a simple treatment of the NL BFKL 
equation, is to emphasize the distinction between the (leading twist) R. G. features, which 
are genuinely perturbative and thus model independent, and the hard Pomeron features 
which will turn out to be strongly dependent on how the effective running coupling is 
smoothed out or cut off at low values of k 2 = 0(A 2 ). 

The problem of the consistency with the R. G. of the BFKL equation with running 
coupling was already analyzed by Collins and Kwiecinski [[5] by introducing a sharp cut- 
off in the transverse momentum integrations. Modifications to the bare Pomeron due 
to the running coupling were also analyzed by Lipatov ||, based on some boundary 
conditions in the soft region which, in our opinion, are eventually equivalent to setting 
a cut-off (see below). The cut-off dependence of BFKL type equations has also been 
investigated |J . Here we whish actually to point out that different ways of smoothing out 
the running coupling in the large distance, small-k region yield different answers for the 
nature, location and strength of the (bare) hard Pomeron, while keeping the validity of 
the leading twist renormalization group factorization, as perhaps to be expected. 

Let us consider the BFKL equation with running coupling introduced in Refs. J7J and 



the corresponding anomalous dimensions ||, 
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where /^(t) = vk^^A(k 2 ) denotes the unintegrated gluon structure function in the hadron 
A, as a function of t = log(k 2 /A 2 ), and the BFKL kernel K(t, t'), which possibly contains 
NL contributions, is supposed to be scale invariant with the spectral representation 



which is also assumed to be symmetrical in t and t', so that X (f + ^0") is even in a. 

The form (1) of the NL equation was proved to be valid for the ^/-dependent part of 
the NL kernel in Ref. ||. But we have also emphasized |TT| that the scale of a s can be 



changed, together with a corresponding change in the scale invariant kernel, so as to leave 
the leading twist solution invariant, at NL level accuracy. Therefore, by assuming Eq. (1) 
we do not emphasize the scale k 2 as the natural scale for the running coupling^, but we 
rather consider a reference form of the BFKL equation which has been widely analyzed 
previously [7-10]. 

We shall also assume that the effective coupling a s (t) = Nc01 ^ & is smoothed out around 
the pole at t — 0, and that the inhomogeneous term fo{t) is peaked at some value t = 
to > 0. We shall often consider in the following the examples 

a 8 (t)~Q(t-i) + ^Q(t-t), (t>0), (3) 

-to), (4) 

but our discussion will not be limited to these particular forms. 

Let us start noticing that if a s (t) < a situ) has a maximum at t = t M , and x{\ > 
x{\) ( as is the case for both the leading and the NL expressions considered so far [[| ||), 
the Pomeron singularity uj p has the upper bound 

uj p < d s (t M )x Qj • ( 5 ) 

This follows from general bounds on the norm of K, and needs no further expla- 
nations. We shall show, however, that besides this general result, the properties of the 

Pomeron singularity are very much dependent on the model for a s (t). 

1 Actually, the NL calculation Q suggests that rather q 2 = (k — k') 2 is the natural scale of the running 
coupling. 



In order to understand this point, we shall discuss the solution to Eq. (1) by using a 
quasi-local approximation of the kernel K valid around the end-point of the 7-spectrum 
in Eq. (2), i.e., 

K(t, f) ~ X Q) (1 + a 2 <9 2 + ....) 6(t - t% a 2 = ±^®. (6) 

This equation is based on the simplest polynomial expansion of xil) around 7 = 1/2, 
which of course has already been used in the literature || . Our purpose here is to obtain, 
by means of the expansion in Eq. (5) a simple picture of the hard Pomeron properties, 
and to show that this picture is stable when some higher order polynomial approximation 
is used. 

By replacing Eq. (6) into Eq. (1), the latter becomes an inhomogeneous second order 
differential equation in the t variable, of the form 

^(/W-/o(t)) = ^(l+a 2 9 t 2 )/(t), (7) 
a s (t) to 

and its homogeneous part is just a Schroedinger-type equation with a given potential V 
and wave number k, given by the expressions 

V(t) =\-(^l-l 



a 2 lu \a s (t) 

^ = (l - |) 4- 0) 

V 00/ a z 

While this potential is always linear for t > i because of the perturbative behaviour of 
the running coupling, its form may vary considerably according to how a s (t) is smoothed 
out in -00 <t <i (Fig. 1 (a), (b)) or cut off (Fig. 1 (c)). 

In the case a s is sharply cut-off (Fig. 1 (c)), one has an infinite potential well, which 
has, of course, a discrete spectrum, whose ground state provides the Pomeron pole, in 
agreement with previous analyses [7-9]. 
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If instead the smoothed out coupling has a flat behaviour below t (Fig. 1 (a)) the 
spectrum is continuum and the Pomeron is a branch cut singularity with branch point 

Up = a s {-oo)x QV (10) 

In this case the eigenfunctions are just plane waves for t — > — oo, and there appears 
to be no reason why their phase-shift should be fixed by some condition M, contrary to 
the case with cut-off. 

Finally in the intermediate case of Fig. 1(b) there may be an isolated singularity too, 
depending on the depth of the well. 

Whatever the relevant case is, the Green's function of the corresponding Schroedinger 
equation is easily calculable, and provides the solution of Eq. (7) if we set f (t) = 
9A$(t — to)- A straightforward analysis shows that for t > t such solution takes the 
factorized form 

V&F A {k 2 ,Q 2 ) = f A (t,u) = f R (t,u)g A t f L {t ,u), (t > t = log ^ (11) 

where /^(/l) denote the regular solution of the homogeneous equation for t — > oo {t — > 
— oo). Their explicit form for t > i is 

f R {t,u) =U{t,u) (12) 
f L (t, w) =f-(t, oo) + R(u)f+(t, uj) (13) 

where R(u) is the reflection coefficient of the well, /+(/_) denote the regular (irregular) 
solutions for t — >• oo in the linear potential, given by the expressions 

W -^MMS?, (t>t), (14) 

Jc± V 2,txi 

x M ./J x( iW = x (i) ( 7 - i) + i x " (i) (7 - i) 3 + .... , (15) 

and C + (C_) denote the regular (irregular) contours for the Airy functions (Fig. 2). 

The reflection coefficient (or 5*-matrix) R(u) in Eq. (13) is easily found - starting from 
k 2 = —y 2 < 0, {u > uj) - for the simple model of a s (t) in Eq. (3). In such case the wave 



functions are just exponentials for t < t, and by the customary matching procedure we 
find 

R(<*) = fr^rA ( 16 ) 

L+{-X ) ~ X 

where L + (L_) denotes the logarithmic derivative at t of / + (/_) in Eq. (14). 

Therefore, the coefficient R(u) contains the Pomeron singularity of the wave number 
X, which in this case is a continuum starting at u — u). It is also clear that R(u) will 
contain the Pomeron singularity for each of the models described in Fig. 1 . In particular 
there may be an isolated pole due to the possible vanishing of the denominator in Eq. 
(16). 

Our final result is thus that for large t, the unintegrated gluon density takes the 
factorized form 

VtfF A (k 2 ,Qf) = f + {u,t)g A t \f-(u,to) + R(u)f + (u,t )) (t > t ). (17) 

The t-dependent factor in Eq. (17), quoted in Eq. (14), is just the nazve regular 
solution |7], 10] of the homogeneous equation corresponding to Eq. (1), for which only 



the perturbative form of a s (t) matters. Its large t, small u behaviour in the anomalous 
dimension regime 

but > X (I) (18) 



2, 

is dominated by a saddle point which reproduces the R. G. behaviour with running cou- 
pling as follows 
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f+(u,t)~ ,^ ex P / lL(a s (t'))dt', (19a) 

vk 2 y/-x J 

1 = ^XU), (19b) 
u 

where jl is the perturbative branch of the anomalous dimension. 

The to-dependent factor consists instead of two terms. One is the naive perturbative 
term f-(u, t ), which is the irregular solution of Eq. (1) and only has, like /+, an essential 



singularity at lu = 0. The other one contains instead the hard Pomeron singularity, due to 
the R(u) factor, which mixes the irregular with the regular solution. It appears, therefore, 
that for large to = log(<5o/A 2 ), the Pomeron term is suppressed with respect to the naive 
term by inverse powers of Ql/A 2 , just because R(u) multiplies the regular solution. This 
explains the mechanism by which the perturbative evolution picture emerges in the case 



of two large scales, i.e. Q > > A , as often assumed in the literature [ff2 



The above transparent picture of the bare hard Pomeron properties emerging from 
the simplest quasi-local approximation of K has the obvious defect of incorporating only 
leading twist solutions in Eq. (19b), due to the use of a quadratic approximation for xil)- 

Let us show that the above picture is only slightly modified when higher twist solutions 
of Eq. (19b) are introduced. We can devise, for instance, a polynomial approximation to 
x{l) °f order 2(n + 1) which has (i) the quadratic part fixed as before and (ii) n higher 
twist solutions of the saddle point equation 

ud-\t) = x in) hn), (20) 

fixed also, together with their mirror solutions 7„ = 1 — 7„, where 70 = 7l(«s(^)) denotes 
the perturbative branch. 

The Green's function of the corresponding 2n + 2-order differential equation will then 
contain n + 1 combinations of type (11), as follows 

n 

Vl^4(k 2 , Ql) = ff(t, u;)g { l>(t , u), (t > t ), (21) 

where the functions f R \t) denote the regular solutions for t — > 00, depending on the 
i th branch of the anomalous dimension, while the functions g^ 1 '(to) are to be determined 
from the continuity (discontinuity) requirements of the Green's function and of its first 
2n derivatives (2n + 1 — th derivative). 

It is not difficult to see that the g^s are the ratio of two determinants of Wronskian 
type involving the f R s and the 's. For large to, they are asymptotically given by 

-1 



gf{to)~ f R l, (t ) , (t G »i). (22) 



r(0, 



6 



We see from Eqs. (21) and (22) that the gluon distribution now contains higher twist 
terms, as expected. Nevertheless, the properties of the leading twist contribution are 
unchanged. In fact, has the behaviour quoted in Eq. (19) and g^ has the decoupling 
property in Eq. (22) for large enough to- Of course the detailed hard Pomeron properties 
determined by the mixing of / + 's and /_'s in the f^s will be more complicated and 
depending on various scattering coefficients. 

The above argument shows that, in general, the structure of the gluon distribution is 
correctly hinted at by the polynomial approximations to x{l)i an d leads us, therefore, to 
the following conclusions. 

Firstly, whatever the weight of the Pomeron is, the solution of Eq. (1) satisfies the 
expected leading twist R. G. factorization property, for any value of to, and precisely 



The t-dependent factor in Eq.(23) is the one predicted by the R.G. , with an additional 
coefficient, which is relevant for NL calculations |5|, pT| . 

Secondly, the hard Pomeron coupled to hadrons is really a nonperturbative phe- 
nomenon, which is very much dependent on the behaviour of the strong coupling in 
the soft region k 2 = 0(A 2 ). Only if the initial scale is large enough does the Pomeron 
decouple, and the BFKL evolution becomes genuinely perturbative. 

If this is the case, then exploring the Pomeron structure and unitarity corrections to 
it becomes really a strong interactions problem in which soft physics plays a major role, 
without much distinction between short and long distance contributions. 

On the other hand, the t-dependent Pomeron singularity present in Eq. (19) will still 
play a role, as a singularity of the anomalous dimension expansion. Therefore, if the 





(23) 



with 




(24) 
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variable a s (t)log(l/x) is not too large, the standard perturbative approach will still be 
applicable, provided a resummation to all orders is performed, along the lines reported 
elsewhere . 
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Figure 1: Form of the potential V ~ a s 1 in Eq. (8), if the behaviour of a s (t) for t < t is 
(a) flat, (b) with a maximum and (c) with a sharp cut-off. 
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